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Abstract 

Starting from the full many body Hamiltonian we derive the leading order energy and density 
asymptotics for the ground state of a dilute, rotating Bose gas in an anharmonic trap in the 'Thomas 
Fermi' (TF) limit when the Gross-Pitaevskii coupling parameter and/or the rotation velocity tend to 
infinity. Although the many-body wave function is expected to have a complicated phase, the leading 
order contribution to the energy can be computed by minimizing a simple functional of the density 
alone. 

1 Introduction 

Rotating Bose-Einstein condensates exhibit fascinating quantum phenomena Uke superfluidity and quan- 
tization of vorticity and their study is currently an active area of both experimental and theoretical 
research. Much of the theoretical work (see, e.g., the monograph [Aj where an extensive list of references 
can be found) is based on an effective description of the ground state in terms of the Gross-Pitaeveskii 
(GP) equation that has recently been proved [LSei IS1| to be exact in a suitable limit. The correspond- 
ing result for the non- rotating case was obtained in |LSYlj . but the rotating case was a long-standing 
problem whose solution required different techniques from those of [LSYlj . In fact, the rotating case 
differs markedly from the non-rotating one since the absolute many-body ground state is in general not 
the same as the bosonic ground state |Slj . 

The limit considered in |LSe[ ISlj is the GP limit of the many-body ground state which means that 
the particle number N tends to infinity while the GP parameter Na, with a the scattering length of the 
interaction potential, as well as the rotational velocity, Q, are kept fixed. In several experiments the GP 
parameter can be quite large, however, and also the rotational velocity can be so large that the effect 
of the rotation becomes comparable with that of the interactions. Some of the interesting phenomena 
expected under such conditions are discussed in P iBPl lECTSO IFl iFBl iFZl [imn limi IHil IWGBPj . To 
describe such cases theoretically it is natural to consider the limit when Na and/or tend to infinity. 
Since the error estimates in jLSe|, IS1| are not uniform in these parameters the results of jLSe| [ST| do 
not apply to this situation and a complementary investigation is called for. In order to allow arbitrarily 
large rotational velocities we shall require that the trap potential increases more rapidly at infinity than 
quadratically (anharmonic traps). 

In a recent paper [CRYl] the limit of GP theory (but not the many-body problem) for large Na 
and large was studied for a rotating 2D gas in a 'fiat' trap with walls that confine the gas to a fixed 
bounded region. (See also |CRY2j for an extension to more general 2D traps.) It was proved that the 
leading order energy and density asymptotics is correctly described by an energy functional of the density 
without a gradient term. This functional for rapidly rotating gases was first introduced and studied in 
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[FBj . Because of its formal analogies with the Thomas-Fermi (TF) theory for fermions it is also referred 
to as a TF functional although the physical situation is quite different. 

In the present paper we derive the TF description from the many-body problem in 3D. This might at 
first sight appear to be a simple combination of the results of [LSei IS1| with those of |CRY1| (extended 
to 3D and more general traps): first take the GP limit of the many-body theory and then the TF limit 
of the GP theory. But such an argument is not valid because the error estimates in |LSe| blow up if 
Na and/or Q tend to infinity. In fact, it is clear that it will at least be necessary to require explicitly 
that the gas is dilute in the sense that the average particle distance is much larger than the scattering 
length. (This condition is automatically fulfilled if Na is kept fixed.) It is also instructive to compare 
with the proof of dimensional reduction of Bose gases in tightly confining traps [LSY41 ISY] where one 
also has the option of taking a limit in two steps. Here the ID or 2D limits of the 3D GP theory do not, 
in fact, cover all cases that can occur when one starts with the 3D many-body theory and takes its tight 
confinement limit directly. In the present situation, however, the TF limit of the GP theory exhausts the 
leading order asymptotics for the many-body ground state, provided only the gas remains dilute as the 
limit is taken. Because the many-body wave function will have a complicated phase it is not self-evident 
that the ground state energy can be computed accurately by a functional of the density alone. That this 
is indeed possible is a basic message of the present paper. 

The main techniques applied in our derivation of the TF limit are on the one hand an extension of 
[CRYl] to 3D and anharmonic, homogeneous potentials, and on the other hand the techniques described in 
[LSY5j for treating the many-body problem in the non-rotating case. Additional tools are the diamagnetic 
inequality as well as the method used in [ST] to bound the many-body energy from above by the GP 
energy in the rotating case. Our results concern the ground state energy and density to leading order in 
small parameters (the reciprocal of the coupling constant or the rotational velocity as well as the ratio 
of a to the mean particle spacing). In contrast to [LSe| Bose Einstein condensation (BEG) is not proved. 
In fact, even in the non-rotating case a proof of BEG in the TF limit is still an open problem. 

The variational wave functions we employ for deriving upper bounds to the energy reflect the expected 
shorts scale structure due to interactions as well as the quantized vortices generated by the rotation. In 
3D the exact vortex structure is presumably rather complicated, e.g., due to bending of vortex lines [Alj . 
Our trial functions do not model such details but they are sufficient for leading order calculations and 
leave room for improvements. 

The paper is organized as follows. In the next section we describe the general setting and state our 
main results about the energy and density asymptotics. Section 3 is concerned with the TF limit of GP 
theory, generalizing the results of [GRYlj to 3D and homogeneous, anharmonic trapping potentials, as 
well as some basic properties of the TF theory needed for the proofs of the quantum mechanical (QM) 
limit theorems. In Section 4 we prove the upper and lower bounds on the many-body quantum-mechanical 
energy and obtain limit theorems for the density as corollaries. Some additional properties of the TF 
density are discussed in an Appendix. 

2 General Setting and Main Results 

We consider a system of N spinless bosons in with mass m, trapped in an external potential V and 
interacting with a nonnegative, radially symmetric pair potential v of finite range. The basic quantum 
mechanical Hamiltonian in a reference frame that rotates with angular velocity Cl is 

N 

i/AT^^ (-A,-4.0 + F(f,))+ (2-1) 

j = l l<i<j<N 

where units have been chosen so that 2m = h = 1. Alsc0 Xj g and Lj = —ixj A Vj for j = 1, ...,N 
are respectively the positions and the angular momentum operators of the particles. The Hamiltonian 
operates on symmetric (bosonic) wave functions in (R'^^, da?i...dxjv) . 

^Wc use the following notation: x = (x, y, z) will always denote a point in K^, while f and z are its projections on the 
X, J/— plane and the 2— axis respectively (cylindrical coordinates), namely x = (r, z). Moreover \x\ and r = |r| will denote the 
moduli of the corresponding vectors, whereas = x + iy will be the complex number associated with the x, j/— coordinates 
of X. 



We keep the external potential fixed and choose the associated length scale (the extension of the 
ground state of —A + V) as a length unit. In order to be able to vary the scattering lengttH of the 



interaction potential v with N we write v{ 



vi{\x\/a) where vi is a potential with scattering 



length 1. Then v has scattering length a and for fixed vi the Hamiltonian is parametrized by N and a 
besides fl. For convenience, we include a factor Att in the definition of the GP parameter 



inNa. 



(2.2) 



The ground state energy of (|2.ip . i.e., the infimum of its spectrum, will be denoted by E'^^ (iV). 

Introducing the vector potential Aq = ■^fl{ez A x) associated with a rotation fl = flcz, where ez is 



the unit vector in z-direction, the Hamiltonian Hf^ can be rewritten in the form 

N 



H 



N 



4" ' ] 



E 

l<i<j<Ar 



V[\X, 



(2.3) 



where r = jez A x\ is the distance from the axis of rotation. The splitting of the term —L ■ VI into the 
contribution of the vector potential and the term — i7^r^/4 corresponds respectively to the Coriolis and 
the centrifugal force in the rotating frame. The vector potential is primarily responsible for the formation 
of vortices while the centrifugal potential affects the overall density profile if the rotational velocity is 
high enough. 

In order that (|2.3p is bounded from below and trapping for all VL we require the external potential V 
to be bounded from below and moreover that 



V{x) 



for 



(2.4) 



Due to the particle repulsion and the centrifugal force the gas cloud expands as 5 or f7 tend to infinity 
so essentially only the behavior of V for large \x\ matters. For simplicity we shall assume that V \s a 
homogeneous functior[f| of order s > 2, i.e. V (Xx) = X'^V (Ax) for any A > and x e M"^, but we need not 
assume that V is symmetric w.r.t. rotations about the z-axis. In order to obtain explicit error estimates 
in Section 3 we shall assume that V is twice continuously differentiable but Holder continuity would in 
fact be sufficient for the main results. 

The case of a 'flat' trapping potential, i.e., V = within a bounded, open set B with a smooth 
boundary and 00 outside, can also be treated by our methods. This case corresponds formally to s = 00 
and the formulas for the energy and density asymptotics can be obtained as the s — > cx) limit of the 
formulas for finite s. At the end of Section 4 we shall comment on the results for flat traps and the minor 
modifications required of the proofs. 

The GP functional in the rotating frame is defined as 



cGP 



dx 



iAc 



V\ 



ln^r'\cl>\'+g\cl,f 



on the domain 



V 



GP _ 



and 



iAn 



e L- 



The corresponding energy is 



inf 



(2.5) 



(2.6) 
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The infimum is, in fact, a minimum and we denote any normalized minimizer by 4>^fi- The corresponding 



density is pf^ 



10^^ p. The minimizer may not be unique because vortices can break rotational 



symmetry, but any minimizer satisfies the variational (GP) equation 



'{V-iAi- 



V 



2ff 



GP^GP 



^Alternatively, one could keep a fixed and vary the external potential and hence the length scale. What matters is only 
the ratio of a to this scale. 

^Asymptotic homogeneity in the sense of |LSY2) would also suffice. 



where is the GP chemical potential. Multiplying (|2.8|) with (/)^^ and integrating gives 

^^';l=E'^J, + 2g\\p^lh. (2.9) 

In [LSet lST] it is proved that E^^{N)/NEfJi ^ 1 as TV -> cx) if .g and Vl are fixed. In the present 
paper we are concerned with the situation where g and/or U, tend to infinity together with N . As we 
shall show, the first term in (|2.5p is negligible in this limit and the ground state energy and density can 
be described in terms of the TF functional 



8jl [p] ^ [ dx{Vp^ \n\^p + gp^] 



(2.10) 



defined on the domain 
with the energy 



P'^P = {p e (r3) : p>Q,Vp(^L^ {M?)] (2.11) 



E''^^ inf £M. (2.12) 



The minimization problem (|2.I2p has a unique solution given by 

pII {£) - Yg ^^"^^^ + ~ ^ + ' ^^-^^^ 

where [•]+ denotes the positive part and mJq is the TF chemical potential determined by the normalization 
IIpJqIIi — 1. Multiplying (|2.I3|) by pJ_Q and integrating gives 

Mro=^.':n+25iiPraiii (2.14) 

By simple rescaling (explained at the beginning of Section 3) we obtain the relations 

^_./(.+3)^TF ^ ^TF ^^^d g-^'^^+'^pll^pll (2-15) 

with 

^^^-(s-2)/(2s+6)^ (2.16) 

and likewise 

9"''^'^pll{9''^'^'^s)^pJl{x). (2.17) 

The case lo — Q corresponds to the standard TF functional without rotation whose relation to the many- 
body problem was already discussed in [LSY5| . Note also that Ej^ is a decreasing function of a; > 
with range (— oo, EjJ^] with E^J^ > 0. In particular, there is an w > such that E];1j = 0. 

In the case that lo tends to infinity the rotational term completely dominates the interaction term. In 
this case it is appropriate to scale lengths with fj^/^""^' rather than (jri/(''+3) (cf. Section 3), obtaining 

^-W(.-2)^TF^^TF ^-2./(.-2)^TF ^^^TF 

with 

^ ^ f^-2(s+3)/(.-2)^ ^ ^-2(.+3)/(s-2) ^2.19) 

and 

^"^'-'^pII {^''^'-'^x) = pZ (f) . (2.20) 
Moreover, as w ^ (X), i.e., 7 ^ 0, we have 

hm i^TF ^ ^TF ^ J f {V{x)-y]<Q (2.21) 

while p^ converges to a measure supported on the set M. of minima of the function W{x) = V{x) — jr^ . 
This together with the other facts mentioned about the TF theory is discussed further in Section 3.2. 

The scaling properties of the TF theory already suggest that one should distinguish between the 
following three cases when the TV — > 00 limit of the many-body ground state with g and/or f2 also 
tending to infinity is considered: 



• Slow or moderate rotation, u <^ 1: The effect of the rotation is negligible to leading order @ 

• Rapid rotation, oj ^ 1: Rotational effects are comparable to those of the interactions. 

• Ultrarapid rotation, uj ^ 1: Rotational effects dominate. 

Moreover, a description in terms of the simple density functional ()2.10p can only be expected in a 
dilute limit. A convenient measure for diluteness turns out to be smallness of the parameter 

a'N\\pYn\\oo-N-'g^pJl\U (2.22) 

Our main result about the energy asymptotics is the following 

Theorem 2.1 (QM energy asymptotics) 

Let V be a homogenous potential of order s > 2, define g = AiraN and lo =^ g^(^^2)/(2s+6)j-^^ ^^j, g^ggn^ne 
that iV~2g3||^TF a.s iV -> oo. 

(i) If g ^ oo and LU ^ as N ^ oo, then lim \ g-'/'^''+^'> N-^E'^^ (N)] = EJ'^ ; 

(ii) If g^ oo and uj>0 is fixed as TV oo, then lim \ g-'^^'+^^ E';^'^ (N)] = Ejl; 
(in) ^ oo and ^ oo as TV ^ oo, then lim \n-^^/^'-^^N-'^E^^ (TV)} = E^f. 

Although stated separately, cases (i) and (ii) can, in fact, be treated together because the convergence of 
the scaled energy is uniform in lo as long as uj stays bounded. 

Besides the energy asymptotics we shall also consider the convergence of the quantum mechanical 
particle density of ground states or approximate ground states. We say that a sequence of bosonic, 
normalized wave functions ^'at G L'^{M.^^) (depending also on the parameters g and fl) is an approximate 
ground state if 

{^iJ^,HN^N)/E^}^{N)^l. (2.23) 

(Recall that E^^ (TV) is, by definition, the infimum of the spectrum of p.ip .) The particle density, 
normalized so that its integral is 1, is defined by 

Pn^q ni^)= / I'^N {x,X2, ■■■,XN)f dx2....dxN- (2.24) 

To ensure convergence of ^ (x) we have to rescale it in accord with (|2.17p . For bounded uj convergence 
of the density follows from the convergence of the energy using standard arguments [Gl ILSi| . 

Theorem 2.2 (QM density asymptotics for uj < oo) 

Under the conditions of Theorem \2.1\ (i) or (ii) we have 

9"^'^''pZ,n (g"''^''^) - pZ i^) (2-25) 

in weak sense. 

The case of ultrarapid rotations, co — » oo, is a little more delicate. Recall that C M'^ was defined as 
the set where the function V{x) — ^r^ attains its minimum. In Lemma 3.2 it is shown that this set is a 
subset of a cylinder and that the scaled TF density ()2.20|1 is eventually concentrated on A1 . In physical 
terms, the ultrastrong centrifugal forces outweigh the repulsion between particles and constrain them to 
the boundary of the available region (in the scaled variables) . We show that the same holds for the scaled 
QM density. 

^This regime could be subdivided further into 'slow' rotations, where vortices do not yet form, and 'moderate' rotations 
where vortices are present. This finer division goes, however, beyond the leading order considerations of the present paper. 



Theorem 2.3 (QM density asymptotics for uj oo) 

Under the conditions of Theorem \2.1\ (Hi) the scaled particle density 

1^*^/(^-2) pQ^ , J (2.26) 

becomes concentrated on the set M as N ^ oo in the sense that the integral of (|2.26|) over any measurable 
set that has a strictly positive distance from M tends to zero. 

The proofs of Theorems I2.1H^751 will be given in several steps. In the next section, we analyze the 
asymptotic properties of the GP and TF theories when g and/or Q tend to infinity. Section 4 contains 
the upper and lower bounds to the quantum mechanical energy that complete the proof of Theorem l2.1l 
Theorem 12.21 is derived from the convergence of the energy by the arguments of [Gl ILSi| and Theorem 
12.31 is proved by showing that any mass outside the set M would be in conflict with Theorem 12. II (iii). 



3 From GP to TF 



In this section we consider the TF limit of GP theory. Since for s < cx) the GP and TF minimizers 
spread out over an increasingly large region as g and/or Q tend to infinity we write x ~ Xx' and 
(j){x) = X^'^/'^cf)' {x') with a suitable length scale A depending on g or so that the relevant x' remain 
essentially bounded. We have ||0'||2 = ll^lb 1, and since V{x) = X^V{x') by assumption, the GP 
functional can be written 



Ax' 



V-iA^.n{S') 0' +X^+^V{x-')\^'\' -\{X^nfr'^\^X (3-1) 



-1 



We now distinguish two cases. When the rotational contribution to the energy is smaller than or at most 
comparable to the interaction energy we equate A^"'"^ with gX^^ , i.e., choose 

X^g^'^'+^\ (3.2) 

For convenience and comparison with [Cmil^myS] and [U we define a small parameter e by 

l/£2 = ^ ^-Ig ^ ^(s+2)/(.+3)^ (3 3) 

The parameter uj = (7~(*~2)/(2'*+^)jl, that measures the relative strength of the rotation with respect to 
the interactions can then be written 

^^^(s-2)/(s+2)f^^ (3.4) 

We now define a rescaled GP functional by writing [(/>] = X^'^E^^ [(/>'], or explicitly, dropping 
the primes. 



[0] 



dx 



V ~iA. 



y\ 



101' + 1-^1' 



(3.5) 



In the next subsection wc study the e — )■ limit of this functional with < oo fixed or tending to zero. 
When the rotation dominates the interaction we take A^^'^ to be equal to (A^fJ)'^ in p.ip . i.e, we take 



Proceeding as before we now write [(/)] = A^^^gp j^/j .^^^.j^ 



V - iAn' 



[v-y]\c^\'+i\ci^t 



where 



^, = ^(.+2)/(s-2) ^^^„2(.+3)/(.-2)^ 

The limit Q oo, lu —^ oo (i.e., Q' — > oo, 7 — > 0) will be considered in Subsection 3.2. 



(3.6) 

(3.7) 
(3.8) 



3.1 The Regime u < oo 

For w < oo fixed or tending to zero, we study the rescaled GP functional introduced in (|3.5p and define 
E?^ by 

The asymptotic behavior of g^^^^^^"^^ E'^fi as g ^ 00 is then given by the hmit of as e — > 0, see 
p.3p . We state this result about E^f^ in the following theorem. 

Theorem 3.1 (GP energy asymptotics for uj < 00) 

Under the conditions of Theorem \2.1\ (i) or (ii) one has, as g ^ 00, 



Proof: The proof is obtained by a comparison between suitable lower and upper bounds for E^^, 
following closely the proof of Theorem 2.1 in [CRYlj . 

The two regimes of slow (i) and rapid (ii) rotations can be treated together. Actually, the lower and 
upper bounds in the case (ii) are sufficient to get the result in the case (i) because the estimates are 
uniform on bounded intervals of m and the TF ground state energy is a continuous function of w. To 
simplify our notation, we denote by C^j any constant independent of e. In particular, Cuj needs not be 
the same from one equation to another, but Ci^ is always uniformly bounded for lo bounded. 

By simply neglecting the positive contribution of the kinetic energy in (|3.5p one obtains the lower 
bound 

e'£?:.[<^]>£ZM']>EZ- (3.11) 
For an upper bound we test the functional (13. 5p with a trial function of the form 



0(f) = Ce^jQeix) Xe{r)ge{r), (3.12) 

where is a phase factor, Xe(r) a function that vanishes at the singularities of g^ and a suitable 
regularization of p^^^ . Note that both and Xe depend only on the 2d coordinate r. For simplicity of 
notation we have suppressed the dependence on lo that is regarded as fixed. 

The regularization of pl^^ is analogous to the one used in [LSY2j . Lemma 2.3 and is explicitly given 
by Qe =je-*^p]:Jj, with 

Je(-)-^exp|-M|. (3.13) 

Since \\je\\i — 1, \/Qs is normalized. It is also clear that converges uniformly to pj^ as e — > and 
it is uniformly bounded in e, i.e., there exists a constant Ci^ such that < C^:- Furthermore, although 
Qe is not compactly supported, it is exponentially small in e for x sufficiently far from the support of 
Pijj- More precisely, denoting 

= sup{|f| : pj^^ix) > 0} < oo, (3.14) 

one has for any a; G M"^, \x\ > 



''UJ I 



We also observe that the gradient of can be bounded in two ways. By using the fact that |Vjj| = £~^je, 
one can easily prove that 

\Vg,\<s-'\g,\ (3.16) 
whereas, by exploiting the regularity of pJJ^, i.e., ||V/9j''j^||^ < C^j, one has 



/ df iV^el < [ dx [ dx' \^pJl{x-x')\j,{x') = ||VpTF|| < 



(3.17) 



^By l|2.13| l the support of pj^ is a compact set, uniformly bounded in uj, if uj remains bounded. 



Writing points r = (x, y) S as complex numbers C^ = x + iy, the phase is defined by 



^(0 - n 



06£ 



C-0 



where £ is a square lattice of spacing defined in the following way: 



L = 



|rj = {m£^,nii,), to, 



r <2R^- 



(3.18) 



(3.19) 



with defined by (|3.14p . We also assume that the spacing is of order y^, i.e., — S-Je, for some 5 > 
independent of e. Note that the phase carries lines of vortices of degree 1 passing through the lattice 
points and it coincides with the trial function chosen in [CRYl] . Moreover the number, N^, of lines of 
vortices included in the support of cj) is bounded by C^/e, because this support is contained in the finite 
region x < R^j. 
The function Xe is given by 



r 1 



Xe{r) = < 



if \r~ rj\ > e'' for all Vj £ £, 



\r - Tj 



(3.20) 



if |r — Tj I < e^, 



for somql 77 > 5/2. 

Finally the constant is fixed by the normalization condition, 110112 = 1, and, since Xe ^ 1 N,, < C/s, 



1 < cl < l + 0(e2')-i) = l + o{e^). 
The rescaled functional evaluated on the trial function p.l2p is given by 



S^^m-ct dx 



dx Qexl 



V - zA, 



/e 9e 



,2 ,0^2 



(3.21) 



(3.22) 



We estimate the three energy contributions separately. Using that (a + 6)^ < 2a^ + 2b^ the first term can 
be bounded bjQ 



dx 



< 2c; / dx 



I df |Vx 



2 



where we have used the trivial bound 



dz gs{x) < C^. 

Denoting by Bl a two-dimensional disc of radius centered at rj E CC] Bu^ , one has 



/ df IVXel' < 
while, by using both p.l6p and (|3.17p . we get 



£2') 



<N,<—, 

£ 



da' 



dx 



IVfe l' - 1 
4& 



< - / df IV^el < — 



Hence p.2ip implies the bound 



df" 



< 



(3.23) 



(3.24) 



(3.25) 



®This requirement on rj is needed in order to apply Theorem 3.1 in |CRY1| (see below in the proof). 
^In this proof Bn, R > 0, will always denote a two-dimensional disc of radius R, centered at the orig 



In order to estimate the second term in p.22p . we first need to restrict the integration to a suitable 
two-dimensional compact set, by exploiting the exponential smallness of Qe given by (j3.15p : 



R3 



\x\<2R^ 



< a 

< a / dr 



dfXe 



132 



\x\>2R^ 

V - iA^/e) 9e 



(v - l^/,) ge 



< 



2 (7 



dx exp 



>2R^ 



< 



V - iA^/e 1 g, 



V - iA^/^ 1 g, 



2 (7 



■ exp 



e 



where 



and we have used the pointwise estimate (see also (3.10) in jCRYlj ). 



which holds for any a? G such that > 2R^. 

The second term of the right hand side of the above expression can be bounded in the following way 
(see also the proof of Theorem 2.1 in [CRYlj ): 



/ 



V - lA^/e ] g. 



< 2 



drxl |V5e|'+2 



A. 



^ - Cnj CuJ Crj , CuJ 



< 



£ e 



4-2r) p3-4j) 



< 



On the other hand we can a ppljEI Theorem 3.1 in [CRYl] to get 



dr 



V - iA^/e 1 g, 



2 . 7ri?2, ^uj 7r\2 C^llogel 



< 



\2 5^ J 



and, by choosing S — ^2-k ju. 



cl / dxQexl 



(V - a^/,) g. 



2 a|loge| 
< H h -5- exp 

e e 



Ru\ ^ a|log£| 



(3.26) 



It remains then to estimate the last term in (I3.22p . namely the TF energy of We have 

and, defining W^(a;) = V{x) — uj^r^/A, 

£Z [Qe] - £Z [pZ] = / * P?,!^ - pII) = / dx p^F (j^ ^W^-W^) 



Differentiabilty of V implies 



\{3e * W^) (x) -WUx)\<^ [ dx' - W^(x)\ e^l-^'l < C^e f dx' r'e~\'^'\ < C^e, 

47r ,/lD.3 ./lD.3 



so that 



(3.27) 



*By scaling the expression can be reduced to an integral over a set contained in a ball of radius 1. The only change with 
respect to Theorem 3.1 in jCRYl) is then the multiplying factor R^. 



Putting all the three estimates (|3.25p . (|3.26p and (|3.27p together we obtain the upper bound 



The upper bound and the lower bound together give the desired result in the rapid rotation regime (ii). 
Indeed, writing as (;-(*-2)/(s+3) g^^^ using (|3.9p . we get, as g ^ oo, 

EZ < g-'^^'^'^EfJ, < EZ + a,g-(^+^)/(2^+«) logg. 

The same estimates prove the convergence in the slow rotation regime (i), by uniformity of the bounds 
and continuity of EJ^. 

□ 

As in the two-dimensional case (see [CRYll ICRY2) V a simple corollary of this result is the L^- 
convergence of the (scaled) GP density \(j)f^\ to the TF minimizer pj^. In particular, in the slow 
rotation regime (i) this gives the convergence of the GP density to pj^^. In that case the proof of the 
upper bound could have been simplified a lot by taking as trial function a suitable regularization of the 
(real) TF minimizer without rotation {uj = 0). 

Also in the rapid rotation regime (ii) other trial functions could be chosen with the same leading order 
contribution to the energy as (|3.12p . One possibility is to cover the support of pj^ with Dirichlet boxes 
and choose a phase factor in each box so that the contribution from A^/^ to the kinetic energy is gauged 
away to leading order. The reason we have chosen a trial function of the form (j3.12p is that the error 
term in p.lOp has the expected dependence on the parameters of a next to leading order term, albeit 
with an unspecified constant in front. In fact, based on considerations of special cases (see, e.g., 0) the 
true GP minimizer is expected to contain, like (|3.12p . a large number ~ uj/e of vortex lines, each 
giving a kinetic contribution of order | logej beyond the leading order term. We note that the optimal 
vortex lattice is expected to be triangular rather than rectangular and bending of lattice lines [Alj will 
also occur but such details only affect the constant factor and higher order contributions in the error 
term in IfJJU^ . 



3.2 The Regime co ^ oo 

We now consider the case of ultrarapid rotations. In the proofs of the quantum mechanical limit theorems 
in Section 4 we shall not make direct use of the GP minimizers to construct trial functions for an 
upper bound to the energy, but the asymptotics of the TF energy and minimizer as cj — > cx3 will be 
important. In this subsection we first derive these properties. The limit theorem for the GP energy for 
ultrarapid rotations can be proved by a simple extension of the corresponding TF result and is included 
for completeness. 

We start with a lemma on the structure of the set M of minimizing points for the sum of the external 
potential V and the centrifugal potential — As always, V is assumed to be homogeneous of order 
s > 2. 

Lemma 3.2 {Ai is a subset of a cylinder) The set A4 of minimizing points for W {x) — y(x) — 
is a compact subset of a cylindrical surface with a fixed radial coordinate 

(3.29) 

where m = — minW > 0. By scaling it follows that all points of the set Mn of minimizing points for 
Wfi{x) = V{x) — i7^r^/4 have the same radial coordinate rn = 0^/'^''~^Vo. 

Proof: That A4 is compact follows from continuity of V and the assumption that W tends to oo as 
\x\ ^ oo. On we have WW = and W = -m, and thus drV = r/2, dzV = and V = r^/4 - to. On 
the other hand, since V is homogeneous of order s, Euler's relation gives rdrV + zdzV — sV and hence 
([3:291) . 



□ 



If V is rotationally symmetric, i.e., V{x) = V{r,z), and strictly monotonously increasing in \z\ 
(examples: V{x) = a|S'|*, or V{x) = ar^ + b\z\'^), then A4 is clearly a circle in the z — plane. It 
is, however also possible that M consists of a two-dimensional subset of the cylinder r = tq (example: 
V{x) = r*/(|z|/r) with / = 1 on some interval but > 1 and increasing outside the interval), or of discrete 
points (example: V{x) = a|x|* + + c\z\'^ with a ^ b). 

Next we consider the convergence of the TF energy E'^^ to Eq^ = inf as 7 = uj^^'^'^+'^yi'^^^^ 0. 

Theorem 3.3 (TF energy and density for a; ^ 00) 

For 7^0 

E^J=EjJ + 0{^y'). (3.30) 
Moreover, the TF minimizer satisfies the hound 

||p^,'il|oo <const.7-3/5 (3.31) 

and for any e > there is a 7^ such that for 7 < 7^ the support of is contained in Ai'^ = {x : 
\x - x'\ < e for all x' G M}. 

Proof: It is clear that E'^\ > Eq^ . For an upper bound we provide a trial function with energy at most 
Eq^ + 0(7^/'^). Let h be any continuous, nonnegative function with support in the unit ball in with 
J h = 1. For S > and a point xq £ M define ps{x) = 6~^h{{x — xo)/S). Then, using that W is and 
that PS is supported in a ball of radius S around xq G Ai, we have 

E^J < S^'^Aps] = I Ax{Wps + ipI]< eJJ + CS' + jS-^h\\l (3.32) 

Choosing S ~ 7^/^ now proves p.30p . 

The TF minimizer is explicitly given by 

Since p'^\ remains normalized as 7 — > 0, it is clear by continuity of W that pj^\ must converge to 
the minimum of W and the support of p!^^ shrinks to the set Ai as stated in the lemma. Moreover, 
E^J > EJJ + 7||p^,^||i and, by (EH, pJJ, = ^JJ, + gWp^lWl and this together with implies 

< [^Z - Wi^)] < l^Z - <i < Oh'/') (3.34) 
and hence ||/c!^_^||oo < const. 7"^^^. 

□ 

Remark: In the case that A4 does not consist of discrete points but is one- or two-dimensional the power 
of 7 in the optimal error terms are of higher order than 7^/^. For instance if V is radially symmetric and 

is a circle, the trial function can be taken to be radially symmetric and the error is 0(7^/^). 

As a complementary statement to Theorem l3.1l we now prove the convergence of the scaled GP ground 
state energy to Eq^ as and a; — s- cx). 

Theorem 3.4 (GP energy asymptotics for uj 00) 

As — > 00 and — )■ 00 

^-2./(.-2)£;GP ^ ^TF ^ Q 1^^,-1 ^ ^2/5^ ^ (3 35) 

with n' ^ and 7 = ^-2(^+3)/(s-2)_ 



Proof: Note first that f]-2V(s-2)^GP ^ n'-'^E^J^ where E^^^ is the ground state energy of the 
scaled GP functional £q^, c.f. p.7p . Dropping the positive kinetic term as in the proof of Theorem 13. II 
immediately gives the lower bound 

^'-'eZ > E^l (3.36) 

For an upper bound we make use of a nonnegative function h with support in the unit ball and J h — 1 
as in the proof of Theorem 13. 3[ this time requiring h to be C°° so that |1Va/7i||2 < oo. For S > we put 
hs{x) = 5~^h{x / 5). Let xq G and define 

0(a;) — \/hs{x — xq) exp{ii7x • {Sz /\xo)/2}. (3.37) 

Testing with this function we obtain 

n'~^E^l<2n'-^V^s\\l+7\\h5\\j+ [ dx(\\ez A ixo-x)\^ + W{x))hs ix~xo). (3.38) 

Since hs{x) = 5~'^h{x/5), the first term is O (fi'^^fJ-^) . Since \\hs\W < ||^5||oo ~ 0((5-3), the second 
term is O (yS~^y In the last integral we use that W ^ C^, that supp/15 is contained in a ball of radius 
6 around xq with W{xo) = E^J, and \ \hs\\i = 1 to get 

/ dx(\ % A (a-o - f)|' + W{x)) hs {x-xo) < E^l + 0{5^). (3.39) 

We thus have 

^^'"'^o,^ < <i + O (r!'-2r 2 + 52 ^ 7^3) . (3.40) 

Equating the second and the last error term leads to the choice b — 7^/^ and an error 
7^/5) — 0(7^/^) provided fi'^^ < 7^/^. For fi'^^ > 7^^^ we choose 8 = il'^^''^ (this corresponds to 
equating the first and the second error term in ()3.40|) ). Then the errors are 0(yt! ). Altogether we 
obtain ([05)1 . 

□ 

Remark: The true GP density is in general not concentrated around a single point in A4 and the trial 
function p.37p is not designed to give optimal error bounds. Another obvious possibility is to replace 
hs by a regularization of the TF density and choose as phase factor of the 'giant vortex' type like in 
|CRY1| . In fact, in the proof of (|4.3p we use a function of this form as an ingredient of the many-body 
trial function. Since Theorem l3.4l is not directly used for the proof of the corresponding many-body result 
we do not elaborate on this point further here. 

4 Proofs of the QM Limit Theorems 

In this section we derive the bounds on the quantum mechanical ground state energy E'^^ that lead to 
the proofs of Theorems I2.1H2.3I The lower bound in the case of ultrarapid rotations is simply obtained 
by dropping positive terms from the Hamiltonian. In the case to < 00 one uses first the diamagnetic 
inequality |LL| to eliminate the vector potential from the Hamiltonian (|2.3p and then proceeds with the 
techniques described in |LSY5j for the non-rotating case. The upper bound for cj < 00 is obtained by first 
bounding the QM energy by the GP energy. The method, that is a generalization of [LSYlj , is described 
briefly in [ST] for fixed g and fi, but in order to keep track of the error terms as g and/or fl tend to 00 
and for completeness we carry it out in more detail. Once a bound in terms of the GP energy has been 
obtained, we can use Theorem 13.11 of the previous section to relate it to the TF energy. In the regime 
of ultrarapid rotation, uj 00, we use a slightly different method that gives an estimate in terms of the 
TF energy and error terms involving directly the TF density whose relevant properties were described in 
Theorem 13.31 The limit Theorems 2.2-2.3 for the density are simple consequences of the energy bounds 
and are discussed in Subsection 4.2. 



4.1 Bounds on the QM energy 

Proposition 4.1 (Lower bound for the QM energy) 

Let the potential V be homogenous of order s > 2. Then 

^-2./(s-2)^-l£;QM ^ ^TF^ (4^^) 

Furthermore, tf uj ^ g-(''-2)/(2s+6)f] fixed and N~'^g^\\pJ^^\\oc ^ as N ^ oo then 

Imiinf {g-^^N-'E^^ (N)} > (4.2) 

uniformly in uj on any bounded interval. 

Proof: To prove (|4.ip consider a normalized A''-particle wave function 5* at and let pn{x) = X'^pn (Ax) 
with A = il^/'*"^ be the corresponding scaled density. Then we can write 

^-2s/(s-2)^-i ^^^^ HN'fN) = inf W, (4.3) 

with W{x) = V (f) - r2/4 and 

C^^ EE r!-2«/(--2) ||[v - an]*^v|[ + J Pn{x) (w{x)- inf dx + 

^j^-2./(s-2)^-l J2 (M'N^vm-X.D^N). (4.4) 

l<i<j<7V 

Since the interaction potential v is by assumption nonnegative the same holds for C.p„ so the left hand 
side of gH) is > inf = E^J. 

Now, let us consider the case when w < oo is fixed as iV — s- oo. By the diamagnetic inequality, 
|(V - iA{x))f(x)\ > |V|/(x)||, cf. [LL], and the bound 

V{x)>p]l + \n\'-2gpYn{£) (4.5) 
that follows from Eq. (|2.13p wc obtain 

Ef^ (N) > Np]J,+ ^ inf^^ Q (*) , (4.6) 

with 

Q(vI/)^^||V.vl/||V ^ v{\x!,-x,\)\^\'d^^x-2gY^ p^n (x.) l^l' d^^ x . (4.7) 

i=l l<i<j<N •' i=l •' 

Since p'^^ix) tends to zero for every x as g — > cxd, cf. ()2.20p it is convenient at this point to carry out a 
rescaling by writing x = Ax' with A ~ g^^^'^+-^\ cf. (|3.2p . The scaled interaction potential w'(x') = A^w(Ax') 
has scattering length a' = y^^/'^'^+'^^a and the corresponding coupling parameter is 

5' = 4^7Va' = g(*+2)/(^+3). (4.8) 

Using (jmi) and ((TT^ we obtain 

g-/(^+3)^-i^QM _ez> f {pzr + wr' ^ inf Q'in (4.9) 

J *, 11*11=1 

where, dropping the primes on the integration variables, 

N „ ^ r 

0'(v|/)^^||V.vI/||V ^ v'i\x,-x,\)\^\'d'^x~2g'Y^ pZ{x!,)\^\'d'^x. (4.10) 

i=l l<i<j<N i=l •' 



We are now exactly in the situation discussed in |LSY5| for the nonrotating case, cf. Eq. (6.61) in [LSY5| . 
Like there, the next step is to divide space into boxes, labeled by a and of side length Z, with Neumann 
boundary conditions and use the lower bound of [LY| for the homogeneous gas in each box. The result 
is (cf. Eq. (6.62) in [LSY5| ) 

g-^^N-'Ef^ (N) -EZ> [ ipZr - E dll'il - Cy""'). (4.11) 

•' a 

Here da is the maximum value of pJw the box a and Y' — a'^N/l^ ^ g'N~'^/P. Since 

9N-'\\pYn\\^=9'N-^pZ\\oo (4.12) 

the diluteness condition implies that Y' ^ for fixed /. If we now first take — > oo and then / ^ 0, 
the Riemann approximation of / (pj^)'^ implies that the right hand side of (|4.1ip tends to zero, proving 
()4.2p for fixed to. It is also clear that all estimates are uniform in uj on bounded sets, so in particular one 
can take uj to 0. 

□ 

Proposition 4.2 (Upper bound on the QM energy for oj < oo) 

Let the potential V be homogenous of order s > 2 and suppose the diluteness condition, -^^^5'^||pJq||oo 
as N ^ oo, is fulfilled. If uj is fixed, then 

hmsup {g-^^N-'E^^ (N)} < E^^ (4.13) 

uniformly in uj on hounded intervals. 

Proof: The proof is a combination of a variational bound on the QM energy in terms of the GP energy 
and the bounds of the GP energy in terms of the TF energy that were discussed in Section 3. For the 
former we can use the same method as in |LSY1] and [S2j and not all details will be repeated here, but 
we shall keep track of the error terms and their dependence on the various parameters. 

The main step is to show that under the stated assumptions 

I, GPii <oil), (4.14) 

by exhibiting a sequence of N particle trial functions ^'^v, iV = 1, 2, . . . , such that 

(^N,H,^r,){^N,^N)-'-NEff, 

AT II GPII —<0{l). (4.15) 

We write the trial functions in the form 

-iiN = F{xi,...,xn)G{xi,...,xn) (4.16) 

with 

G{x,,...,XN)^nll^fxli^^) (4.17) 
and a real function F. Partial integration, using the variational equation (I2.8P and the reality of F, leads 
to 

J2 I |V.Fp|Gp+ I v{\x,-x,\)\F\''\G\^-2g ^ / p'^lmF?\G? ■ (4.18) 

The second line of (|4.18p is a real quadratic form in F and we shall make use of the fact for an upper bound 
on the bosonic ground state energy is not necessary to require that the trial function F is symmetric 



under permutations of the variables. This can be seen by a simple adaption of an argument of Lieb [Llj 
which implies that the infimum over all functions F is the same as the infimum over all nonnegative, 
symmetric functions. 

Like in [LSYlj we shall take a trial function of the Dyson form [D] 

F{xi, . . . ,xn) = UZlF^i^u ■ ■ ■ ^^i) (4-19) 

where 

Fi{xi,...,Xi) = f{ti), ti = min(|xi - Xjl, j = 1, . . . ,z - 1) , (4.20) 
with a function / satisfying 

< / < 1, /' > 0. 
The function / will be specified shortly. Our estimates involve the quantities 

f (l-f), J^l (f + \vf), ff. (4.21) 

By exactly the same computation as leads to Eq. (3.29) in [LSYlj we obtain, provided -/V||/3?olloo^ < 1' 



I JR3« 7r3JV 



' 't2 T f „GP/-f\2 I 2 3 2|i GP||2 



and the same technique gives also a bound on the last term in ()4.18p 

-2g\\FG\\^' I Ps':o(^0m'|Gp<-2giv/ dxp^l{xf + 2gN^I\\p%\\1. (4.23) 

We now choose the function /. For a parameter h > a that will soon be fixed we define 

/(.) = /(l + ^i)"WA forr<& ^^24) 
1 1 for r > b 

where u{r) is the solution of the scattering equation 

-u"{r) + \v{r)u{r) = with m(0) = 0, lim u' [r) = 1 

r — >oo 

and ei is determined by requiring / to be continuous. Convexity of u gives 

I for r < a , , , I for r < a 

r > u{r) >\ - , 1 > u'ir) >{ " 

r — a for r > a 1 — - lor r > a . 



These estimates imply 



J<(l + ei)47ra (4.25) 

/ < 47r + afe(6 - a)^ (4.26) 

K < 4:7r{l + ei)a (^b ~ (4.27) 

< ei < — (4.28) 
b — a 



Before proceeding further, we need to relate the supremum of the GP density, Hp^oIIoc IIpJ,oIIoo 
since the diluteness condition is stated in terms of the latter. For this purpose we write = Re^^ with 
real S and the nonnegative amplitude R. A straightforward computation, using V • A = 0, gives 

= -{AR)e"^ - 2i{VR) ■ {VS)e'^ + {VSfRe'^ - i{AS)Re''^ 
+2iA- (Vi?)e''^ - 2A- ("^ S)Re'^ + A^Re'^ 



and from the GP equation (12. 8|) one obtains 

-A + [VSf 2A- (V5) +A^ + 2.gpGP + y _ if^2^2^ ^ ^ ^gp^_ 

In any point x e M'^ where R is maximal p^^{x) = ||p^j^||oo and Ai?(x) < 0. Thus, 

2g||p^Fjoo < '{^S{x)f + 2i*(a;) • V^(x) - A^x) - V{x) + {^fr^ + p^P . 
and since -{VSf + 2A-VS < A^ we obtain 

25l|p^l^||oo < - mf^ {V{x) - . 

On the other hand, by (|2?T3l) we have 

25l|pJ.!;iloo = pj,?,- mf {V (x) - (4.29) 

and therefore, using (|2.9|) and (|2.14p . 

H\P%\U < M^o - pJ,o + 25l|pJ,nlU < - EjJ, + g\\p%\\^ + 2ff||pJ,f,|U, (4.30) 

i.e. 

1 + ■ (4.31) 

By using (|2.15p and Theorem 13. II one sees that 

II TP II - llPl.t^jlloo (ff Eg,n~ El,uj - 0(1), (4.32) 

so 

\\p%\\^<2\\p^l\\oo{l + o{l)}. (4.33) 

The diluteness condition IpJqI |oo thus implies the corresponding condition for the GP density, 

i.e., ^~^ff'^||p^nlloo — > 0. Therefore, by choosing 

b=impfJi\U-K (4.34) 

it follows from K^ - KM that ei ^ 0, £3 = A^IIp^FJccJ" -> and eg = g-^N^K^\\p'^^\\oc for 
^ cx). Altogether one gets from (gTH), using ((T^ . ' (li!^ and that 

Ar-ii?Q^ < i?,°r^ + <?llp,^Fll2{0(£i) + 0(e2)} + 5llp,^FllooO(e3) < + o(l)<?||p«P |U, 
i.e. (mm). By ((TTT)) we have 

9''/^'^'^ {Sllprolloo} = ff^/(^+^)||pJ,!;iU = IIpI^IIoo, (4.35) 
and by (|4.33p and Theorem 13. II we can conclude that 

3-/(^+3)^-l^QM < 5-./(.+3)^GP + ,(1)||^TP ^ ^TF ^ ^ (4 gg) 



□ 



For ultrarapid rotations the proof of (|4.33p given above is not valid because the error term may blow up 
as w ^ oo. We shall therefore treat this case separately, using a trial function different from (j4.17|l . If 
a general proof of ||p^olloo ^ (const.) ||pJq||cx) can be found, then Eq. ()4.15|) is verified also for w ^ oo 
and the proof of the next proposition would follow in the same way as the previous one. 

Proposition 4.3 (Upper bound on the QM energy for lo oo) 

Let the potential V be homogenous of order s > 2 and suppose the diluteness condition, N^'^g^\\p^\f\\oo 
as N —^ oo, is fulfilled. If Q —^ oo and cj — > oo as N oo, then 

limsup {n-'^/(^-'^N-^Er^^{N)} < EJI 



Proof: The first step is to choose a suitable phase factor for the trial function to compensate the vector 
potential in the kinetic term as far as possible. As shown in Lemma 3.2 the set of minimizers of 



Wn {x) = V (x) - 



is a subset of a cylinder with radius ro > 0. We define the phase factor as follows: 

N 



'd{xi,...,XN) = ]^W(a;i 



with 



(X ) = exp < I 



i9 > for af = (r, i?, z), 



(4.37) 



(4.38) 



(4.39) 



and where [•] stands for the integer part. Any function 'i' E can be written as ^' = G$ with $ e 
and a straightforward computation gives 



N,g,n 



(*) 



- «o($)+E/ 
2VRe(/ $ 



2J2. 



n'r 



l<i>l 



-iVj-Ao(f,) e |$| 



((-*v,-Ao(x,)) e) v,a>' 



Since 



(- 



iV - A< 









rn\ 






2 " 





we get, using the Cauchy- Schwartz inequality combined with 2ab < + 



N 



I<i>l 



2 " 



rjQ 
2 



1*1 



|V,$| 



In particular, since [^^'^^^j = + with |k| < 1, we then have 



'QM 



with 



'QM 
N,gS1 



N 



2|V,$|- 



E 

l<i<j<N 



vi\x.-x,m' 



(4.40) 



(4.41) 



(4.42) 



(4.43) 



(4.44) 



and 



The functional (|4.44p describes the QM energy of a non-rotating system of particles with mass 1/2 in the 
positive external potential {Wq — inf M^n) and with a two-body interaction potential v >0. 

We need to chose a trial function for (|4.44l) . As as in the proof of the previous proposition we are 
dealing with a real quadratic form so the infimum over all <i> G L^(M'') is the same as the infimum over 
all symmetric $ e L'^{R^). Thus, by 

EZ^-NMWn<M'-^^0^. (4.46) 

We denote by pe = je * pj^n the regularized TF density with defined by p.l3p for any e > 00 Observe 
that ll/Oflli = 1 and ||Vy^||2 = since |Vje| — e^^je and = 1 for any e > 0. Our trial function 

is defined as 

...,xn) = F{xi, ...,fjv)G(fi, ...,xn), (4.47) 

with 

N 



G{xi,...,XN)=Y[^Peixj) (4.48) 

while the function F is of the Dyson form, cf. ([4T9| - (|420|l and (|4?24l) . 

The estimation of f ;^|^fj(l>)/| |l>| |^ follows closely the computations in [LSYlj . Eqs. (3.11)-(3.29), but 
with the regularized TF density pe instead of the GP density. The diluteness condition ^^^5'^||pJ^ol|oo 
implies that the same condition is also fulfilled for p^ because 



\\Pe\U<\\p','n\\oo\Mi = \\p',!n\U (4.49) 

The GP equation, that was used in the computations in [LSYlj to obtain Eq. (3.28) in that paper, is not 
at our disposal for but we can instead use the Cauchy-Schwarz inequality combined with 2ab < + b^. 
In this way, using also the positivity of w, we obtain 

ETgn-NinfWn < ^ 4 / \V,G\' + (Wn{x,)^miWn)F'G 

Since = ^~^je we have for the first term of (|4.50p the bound 

4||^^G||-' / \V,G\^F^ = A\\FG\\f [ 



v{\x,-Xj\)F^G^} . (4.50) 

R3" 



p, (Xj) 



VjPe (Xj) 



2y^p, (Xj) 



F^G^ < \. (4.51) 



Using Eqs. (3.15)-(3.16) and (3.21) in |LSYlj as weU as (Wn - MWn) > 0, we obtain the bound 



\\FG\ 



-2 



/ (Wn (xj) - inf Wn) F^G^ < }. ,, ^ / d^x: (Wn - inf Wn) Pe (4.52) 

V / i — iV ||Pe||oo-f Jr3 V / 



with / defined by (|4.2ip . provided 7V| joo-/^ < 1 that is guaranteed by the diluteness condition. The 
last two terms in (|4.50p are bounded in exactly the same way which leads to (|4.22p with in the place 
of pffi- We omit the details. Altogether we have the upper bound 

N-^E'^^^-MWn<^ + il + o{l)) f dj:{(Wn~infWn)pe + 4.gp'i}+oil)g\\pe\\oo- (4.53) 
Since pe is normalized, 

||Pe||2 < IIp.IIoo < ||pJ,?;i|oo (4.54) 



^This e is unrelated to the e defined in 113. 3| l. 



by (|4^ . Also, 

n-^^/i^-^)MWn^EJl (4.55) 

..... . . 2 

and by Lemma 13. 21 all minimizing points have the same radial coordinate rji — rQfl^-^ with rp the radius 
of the set M . The inequality (|4.53p now implies 

+5i7-2^/(^-2)^iipT?;iu(i+o(i)) 

^^-2s/(.-2)g-2^ (4.56) 

We now have to bound each term of the right hand side. For the second term we use that 

n-''/^'-'^g\\pj:'n\\oo=7\\pZ\\oo <oil) (4.57) 
by (|2.20p and because ||p!^^||oo < 0(7"'^/^) by Lemma 3.2. The last term in (|4.56p is o(l) as long as 

It remains to consider the first term in (|4.56p . By scaling we have 

Q-2s/is-2) f _ jj^f ^\ ^ f dxlv (x) - y + 

(4.58) 

with 

Pe(x) = n^p, [n^x^ = jr*p^F(f), i=n-Ae. (4.59) 

As 7 ^ the support of p]^^ becomes concentrated on the set A4 where V (x) — jr^ is minimized and all 
points in M have the same radial coordinate, rp. Moreover, outside of the support of pl^^ the regularized 
density p^ decreases exponentially if d{x) /? — > cxd where d{x) is the distance of x from the support of 
p'J^\ (see also the proof of Theorem 3.1). This implies that the right hand side of (|4.58p tends to zero as 
7 ^ 0, O ^ oo and e ^ with e-if]-'*/(''-2) o. 

□ 



f^/2^2 



pTF 



4.2 Convergence of the QM particle density 
4.2.1 The case uj < oo 

Proof of Theorem \2.S[ We use Griffiths' argument [G] in the same way as for an analogous problem in 
[LSi| . Take any bounded function f : M.^ ^ M. and for any a € [S, 5] {S > 0) perturb the Hamiltonian 
Hpf with the external potential ag'^^'^^'^^^ f{g~^^^''~^^^x). Because / is bounded, the statements (i)-(ii) of 
Theorem 12.11 can also be proven for the perturbed external potential {V (x) + (Tg'*/^^+"^'/(g^^/('*+"^'a;)}. 
Namely, if g — > oo and w > is fixed then the corresponding ground state e'^^ ^ (N) converges to 



N 



Um^ {g-'^^'+'^N-'E^^^^ (N)} = EJl^, for any a G [5, <5] , <5 > (4.60) 



where Ejl „ is the TF energy with V replaced by + af. Consequently, for any approximated ground 
states \1/Ar of Hm we have 

(*jv, [Hn^. - Hn] ^n) = Ng'/'-'+^^a (/, pn) , (4.61) 
with pn{x) = (7'^/(*+^^pAr(g^/('^+'^)x). Hence, the Rayleigh-Ritz principle for any a € [S,S] leads to 

Ef^,. (N) - i?,^o.o W < ^ff^/'^+^V {f,pN) + {^N, [Hn E^^ (N)] vf^) . (4.62) 



Because (fSJS]) and (|2T7| are still valid with V{x) replaced by V{x)+ag''^''''+^^ f{g-'^^^^+^^x), the previous 
inequality combined with (j4.60p implies in the limit N oo that 



<^'^ [EZ.a - EZ.a] <liminf (/, p^) for any a € (0, S] , (4.63) 

■ ■ ■ ■ A' — > + oo 

whereas a negative parameter a G [—(5, 0) reverses the inequality. The proof of the differentiability of 
at (t = is deduced from similar estimations as from (|4.6ip to ()4.63p combined with the continuity 
of the function a i—^ (^f, pj^ , where pj^ ^ is the minimizer of the variational problem Ej^ ^. We 
omit the details. In other words, we have dcEj^^^ — {f,pjjj). Consequently, by (I4.63|) and its reversed 
inequality, we obtain Theorem 12.21 

□ 

4.2.2 The case — > oo 

In the case of ultrarapid rotations Griffiths' argument is not as easily applicable as in the previous 
situation. There are two complications. First, perturbing V with a scaled additional term af leads to 
a potential that is not homogeneous and the proof of the upper bound for the QM energy has to be 
modified in order to get (|4.63p with Ej\ ^ replacing Ej^ ^. Secondly, if M. consists of more than one 
point the variational problem for Eq\ ^ does not have a unique minimizer and Eq\ ^ is not differentiable 
at tr = in general. To get around this complication, we would need to perturb the Hamiltonian with 
an additional term to avoid the degeneracy of the variational problem. An example of such an argument 
is given in [LSe] . But since we are content with proving that the density is concentrated on M. in the 
limit and not striving to obtain the exact limiting measure on AA we can ignore all these problems and 
use the following simpler argument. 

Proof of Theorem \2.S[ If K, is any set with a positive distance from A4, then W{x) — Eq^ > c > on 
JC with some strictly positive number c. Hence, by Eqs. (|4.3p - (|4.4p . 

^-2s/(.-2)^-i ^q,^^HN'fN) - E^^, >c I pN. (4.64) 

JK 

On the other hand, by Proposition 4.3 we know that the left hand side of (|4.64p tends to zero as iV — > oo, 
so pn 0. 

□ 

4.3 Remarks on 'flat' trapping potentials 

As mentioned in the Section 2 the case of a 'flat' trap, that corresponds formally to s = oo, can be treated 
in essentially the same way as we have done for s < oo. By a flat potential we mean that is inside 
some open, bounded set B with a smooth boundary and oo outside. More precisely, the kinetic term in 
the many-body Hamiltonian (12. 3p and the GP functional (|2.5p are defined with Dirichlet conditions on 
the boundary of but Neumann conditions lead, in fact, to the same results in the large g, large 17 limit. 

For s = oo Eq. (|2.16p resp. (j2.19p reduces to cj = g~^/'^Q. resp. 7 = f7~^.g and the TF energy scales as 
g-'E-^l = Ej%, resp. ^'^E^J, = E^^. Theorem 2.1 becomes in case (i) {g-'N-^E'^^ {N)] ^ i^^^, 

in case (ii) [g-'N-^E^^^ (iV)} ^ Ej^, and in case (iii) [n-^N'^E'^^ {N)] - E-^J. 

The proofs require only some minor modifications. For instance, in the proofs of Therorem 3.1 and 
Propositions 4.2-4.3 the trial functions have to be modified in order to take the boundary condition into 
account. In the case of ultrarapid rotations it has to be noted that the set M is now a subset of the 
boundary dB of B, consisting of the points on dB where the centrifugal potential — r^/4 is minimal, i.e., 
where r is maximal. If the boundary is smooth one can still use a Taylor expansion for the proof of an 
analogue of Theorem 3.3 (with different error terms). 



Appendix A The TF density at large rotational velocities 

The TF density is explicitly given by (|2.13|) . i.e., 

pjl i^) = ^ [f^Yn + - y i^)] + • (A.l) 

To get a picture how it changes with the parameters, in particular as uj = (^^(''^2)/(2s+6) _^ jg 
convenient to use the scaling 

17«/(-^)pJJ, (n'/is-^^x) = pZ ix) (A.2) 
to eliminate the dependence of the potential on il and consider 

with 7 = u;~2(s+3)/(s-2)_ ^ incrcascs from to oo, 7 decreases from +00 to zero and, due to the nor- 
malization, the chemical potential ii^^ decreases monotonically from +00 to Eq^ = min {V (a?) — jt"^}- 
Since V is homogeneous of order s > 2 it is clear that Eq^ < 0. By continuity there is a jc (and a 
corresponding ojc) such that fJ^l^fi = 0. Explicitly, 

^, = lJ[y-V{x)]^dx. (A.4) 

Since V(0) = and V is continuous, we have Pj^{x) > in a neighborhood of for 7 > 7c (i.e., uj < uJc)- 
For 7 < 7c, on the other hand, pj^^i{x) = in a cylinder around the z axis, because /i^^i < and V > 0. 
In other words, for 7 < 7c (i.e., ui > ujc), the centrifugal force creates a 'hole' of radius > 2-^/^^/1x1 in the 
density. 

To describe this in a little more detail we introduce cylindrical coordinates (r, z, z9) and consider the 
density as a function of r at fixed (z, 1)), first for z = 0. Using homogeneity of V we see that the boundary 
of the support of p!^^ is determined by the solutions to the equation 

a{^y - rV4 = (A.5) 

with a(i?) — V{1, 0, -d) > 0. If 7 > 7c (i.e., uj < lUc) so pl^^i > 0, the equation has one solution, r(i9)+ > 0, 
and pj^i{r, 0, -i?) > for r < r(i?)+ but pl^^i = outside. Eq. (jA.5|) has two solutions, r{{))^, if 7c > 7 > 7^ 
(i.e., LOc < < uj^) where 7^ > is determined by solving (|A.5P together with 

sa{^y~^ (A.6) 

The density p^_^(r,0,?9) is nonzero for r inside the interval (r('i9)^ , r(z?)+) but vanishes outside. For 
7 = 7ij the interval shrinks to a point with radial coordinate 

r(i?)'"" = (2sa(i9))-i/("-2) ^ gyp ^ ij^f ^(^)+ (A 7) 

7>7a 7>70 

The chemical potential ijJ^Ji and hence 7^ is determined by inserting (|A.7|) into (|A.5|) . If 7,3 > 0, then 
p'^\ is identically zero in the t9 direction for 7^ > 7 > 0. If 7,3 = 0, on the other hand, then (jA.5|) implies 
that {rf'^,Q,'d) E M and thus r(i?)^™ = vq by LemmaEH 

For z ^ we can use the homogeneity of V to write V{x) = r''V{l,z/r, d)a and apply the considera- 
tions above to a ray with fixed z/r and 'd. It is, however, more interesting to consider what happens at 
fixed z 7^ in the case that V is monotonically increasing in |z|, for instance if V{x) = Vo(r) + c\z\'^ . At 
fixed z the term c|z|* acts as a shift of the chemical potential, pj^^^ — c\z\^. Hence for a 'hole' to 

appear at z > it is not necessary that becomes negative, it appears already when p']^^ — clz]'* < 0. 
For any 7 < 00 the density has therefore a 'hole' for sufficiently large |z|, the width increasing with |z|. 
As 7 decreases from an initial value > 7c the 'hole' moves down to lower values of z, reaching z = at 
7 = 7c- 



As an example, eonsider V{r, z, t9) = r*(l— e sin^ + with < £ < 1. Here 7^1 > for all ?9 G [0, 2tt) 
except for i9 = and = tt. As 7 ^ the density converges to {^(5(a; — 1) + ^^{x + 1)} S{y)S{z). If 
e = 0, i.e., V = + c\z^\, then 7,? = for all and the limiting density for 7 — > is (2^) 6{r — l)6{z). 

Acknowledgments. This work was supported by an Austrian Science Fund (FWF) grant P17176-N02. 
JY would like to thank the Niels Bohr International Academy and Nordita, Copenhagen, for hospitality 
and Chris Pethick, Gentaro Watanabe and Gordon Baym for discussions. 



References 

[A] A. Aftalion, Vortices in Bose-Einstein Condensates, Birkhauser, 2006. 

[Al] A. Aftalion, T. Riviere Vortex energy and vortex bending for a rotating Bose-Einstein conden- 
sate, Phys. Rev. A 64 (2001), 043611. 

[B] G. Baym, Rapidly Rotating Bose-Einstein Condensates, J. Low Temp. Phys. 138 (2005), 601-610. 

[BP] G. Baym, C.J. Pethick, Vortex Core Structure and Global Properties of Rapidly Rotating Bose- 
Einstein Condensates, Phys. Rev. A 69 (2004), 043619. 

[CRYl] M. CORREGGi, T. Rindler-Daller, J. Yngvason, Rapidly Rotating Bose-Einstein Conden- 
sates in Strongly Anharmonic Traps, ArXiv math-ph/0606058, J. Math. Phys. (in press). 

[CRY2] M. Correggi, T. Rindler-Daller, J. Yngvason, Rapidly Rotating Bose-Einstein Conden- 
sates in Homogeneous Traps. Preprint 

[D] F. Dyson, Ground-State Energy of a Hard Sphere Gas, Phys. Rev. (1957), 20-26. 

[ECHSC] P. Engels, I. CoDDiNGTON, P.C. Haljan, v. Schweikhardt, E.A. Cornell, Observa- 
tion of Long-lived Vortex Aggregates in Rapidly Rotating Bose-Einstein Condensates, Phys. Rev. 
Lett. 90 (2003), 170405. 

[F] A.L. Fetter, Rotating Vortex Lattice in a Bose-Einstein Condensate Trapped in Combined 
Quadratic and Quartic Radial Potentials, Phys. Rev. A 64 (2001), 063608. 

[FB] U.R. Fischer, G. Baym, Vortex States of Rapidly Rotating Dilute Bose-Einstein Condensates, 
Phys. Rev. Lett. 90 (2003), 140402. 

[FZ] H. Fu, E. Zaremba, Transition to the Giant Vortex State in a Harmonic-plus-Quartic Trap, Phys. 
Rev. A 73 (2006), 013614. 

[G] R. Griffiths, A Proof that the Free Energy of a Spin System is Extensive, J. Math. Phys. 5 (1964), 
1215. 

[KTU] K. Kasamatsu, M. Tsubota, M. Ueda, Giant Hole and Circular Superflow in a Fast Rotating 
Bose-Einstein Condensate, Phys. Rev. A 66 (2002), 053606. 

[KB] G.M. Kavoulakis, G. Baym, Rapidly Rotating Bose-Einstein Condensates in Anharmonic Po- 
tentials, New J. Phys. 5 (2003), 51. 

[LI] E.H. LiEB, The Stability of Matter: From Atoms to Stars, Bull. Am. Math. Soc. 22 (1990), 1-49. 

[L2] E.H. LiEB, Thomas-fermi and related theories of atoms and molecules, Rev. Mod. Phys. 53 (1981), 
603-641. 

[LL] E.H. LiEB, M. Loss, Analysis (2nd. ed.), American Mathematical Society (2001). 

[LSe] E.H. LiEB, R. Seiringer, Derivation of the Gross-Pitaevskii Equation for Rotating Bose Gases, 
Commun. Math. Phys. 264 (2006), 505-537. 



[LSi] E.H. LiEB, B. Simon, The Thomas-Fermi theory of atoms, molecules and solids, Adv. Math. 23 
(1977), 22-116. 

[LSYl] E.H. LiEB, R. Seiringer, J. Yngvason, Bosons in a Trap: A Rigorous Derivation of the 
Gross-Pitaevskii Energy Functional, Phys. Rev. A 61 (2000), 0436021-13. 

[LSY2] E.H. LiEB, R. Seiringer, J. Yngvason, A Rigorous Derivation of the Gross-Pitaevskii Energy 
Functional for a Two-dimensional Bose Gas, Commun. Math. Phys. 224 (2001), 17-31. 

[LSY3] E.H. LiEB, R. Seiringer, J. Yngvason, One-Dimensional Bosons in Three-Dimensional Traps, 
Phys. Rev. Lett. 91 (2003), 1504011-4. 

[LSY4] E.H. LiEB, R. Seiringer, J. Yngvason, One-Dimensional Behavior of Dilute, Trapped Bose 

Gases, Commun. Math. Phys. 244 (2004), 347-393. 

[LSY5] E.H. LiEB, R. Seiringer, J. P. Solovej, J. Yngvason, The Mathematics of the Bose Gas 
and its Condensation, Oberwolfach Seminars 34, Birkhauser (2005). 

[LY] E.H. LiEB, J. Yngvason, Ground State Energy of the Low Density Bose Gas, Phys. Rev. Lett. 80 
(1998), 2504-2507. 

[Lu] E. Lundh, Multiply Quantized Vortices in Trapped Bose-Einstein Condensates, Phys. Rev. A 65 

(2002) , 043604. 

[SY] K. SCHNEE, J. Yngvason, Bosons in Disc-Shaped Traps: From 3D to 2D, Commun. Math. Phys. 
269 (2007), 659-691. 

[SI] R. Seiringer, Ground State Asymptotics of a Dilute, Rotating Gas, J. Phys. A: Math. Gen. 36 

(2003) , 9755-9778. 

[S2] R. Seiringer, Dilute, Trapped Bose Gases and Bose-Einstein Condensation, in Large Coulomb 
Systems, Lect. Notes Phys. 695, 251-276, J. Derezinski, H. Siedentop, eds.. Springer (2006). 

[WGBP] G. Watanabe, G. Baym, S.A. Gifford, C.J. Pethick, Global Structure of vortices in 
rotating Bose-Einstein condensates, Phys. Rev. A, 74 (2006), 063621. 



